Virtual photon processes in the collective dynamics of multilevel atomic ensembles 
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We study the contributions of virtual photon processes to the dynamics of a system of A'^ multilevel 
atoms sharing one excitation and interacting with the quantized vector electromagnetic field. The 
Rotating Wave Approximation significantly simplifies the derivation of the equations of motion 
describing the collective atomic dynamics, but it leads to an incorrect expression for the dispersive 
part of the atom-atom interaction terms. For the case of two-level atoms and a scalar electromagnetic 
field, it turns out that the atom-atom interaction can be recovered correctly if integrals over the 
photon mode frequencies are extended to incorporate negative values. We explicitly derive the 
atom-atom interaction for multi-level atoms, coupled to the full vector electromagnetic field, and 
we recover also in this general case the validity of the results obtained by the extension to negative 
frequencies of the formulas derived with the Rotating Wave Approximation. 



I. INTRODUCTION 

Since the early work of Dicke [l| , the problem of collec- 
tive emission from a sample of identical two-level atoms 
has received much attention, see e.g. 0,01 and references 
therein. The problem is typically treated using the Ro- 
tating Wave Approximation (RWA) to mathematically 
simplify the derivation. A recent series of papers [3t3] 
has discussed the validity of RWA and, particularly, the 
importance of effects due to virtual photons in a system 
of two-level atoms and a scalar light field. The equations 
of motion without RWA for a vector field interacting with 
a system of two-level atoms were derived by Milonni et al. 
0, who identified differences between the RWA and the 
full equations of motion. They also found that the atom- 
atom interaction terms which were not obtained correctly 
by the RWA derivation, can be obtained if integrals over 
photon mode frequencies in the RWA expression are ex- 
tended to negative values. This artificial extension of 
the integration range leads to the inclusion of rapidly 
oscillating terms which, indeed, contrast the argument 
for applying the RWA in the first place. The artificial, 
negative energy photon states, however, formally violate 
energy conservation by the same amount as the far-off 
resonant, virtual states, which may really influence the 
atomic evolution. This may explain why they lead to the 
same dispersive coupling among the atoms. 

The equivalence of the two calculations has been 
demonstrated by going through both calculations and 
comparing the final results for the two-level atoms with 
scalar fields. It has not however been proven by a physical 
argument. Hence it is not clear if a similar simplification 
holds in the general case of multi-level atoms and vec- 
tor electromagnetic fields. Although the full equations of 
motion without the RWA in this general case have been 
derived by Friedberg et al. [IQ], the question about the 
equivalence of this result to the extended RWA derivation 
remains open. In this work we explicitly derive the equa- 
tions following from the usual RWA and from the RWA 
with the extension to negative frequencies and present 



the parallel to the full equations. Our results confirm, 
also for this general case, the validity of the much sim- 
pler extended RWA ansatz. 

The paper is organized as follows. In Sec. II we derive 
the general multi-atom collective dynamics within the 
RWA and within the RWA extended to negative frequen- 
cies. In Sec. Ill we derive the general multi-atom dy- 
namics, including virtual excitation of non-RWA dipole 
coupled states with one photon and two excited atoms. 
Sec. IV summarizes and concludes the paper. Evalua- 
tion of a few matrix elements and integrals are provided 
in two appendices. 



II. ROTATING WAVE APPROXIMATION 

We consider a system of N atoms located at rest at 
positions Vj. Each atom has three degenerate excited 
states |e^"'^), je*^) and |e+^) with excitation energy hujQ 
above a single ground state \g) and a dipole moment deg- 
The three excited states decay with the same rate F = 
4<ieg^ and emit photons with the polarization dictated 



by dipole selection rules 

The Hamiltonian describini 
Schrodinger picture is [1, [HI, [I 



the system in the 
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FIG. 1: (Color online) Couplings due to the interaction 
Hamiltonian Hv- The system with initially only one atomic 
excitation is directly coupled to the atomic ground state with 
one photon emitted and virtually coupled to a state with two 
atomic excitations and one photon. In the RWA the virtual 
coupling is neglected. 



Taking the transition dipole moment operator for atom j 
D'=deMg+cflJ (4) 

and the electric field operator evaluated at the position 
of this atom 



V J 



(5) 



we can write the interaction Hamiltonian in the form 

N 

i=i fc,A 



^k\^k\^ 



(6) 



Here at^ denotes the photon creation operator in a 

mode with wave vector k, energy lj^ — ck and po- 
larization along direction ej:^, A — 1,2. The vacuum 

d,,n(^)K The dipole op- 



coupling constant is gk = deg (^f^) 
erators for atom j are dig 



J2^=-idug\ej){gj\ and 
^ge — J2l=-i'^gi^\9j){^'j\ with the complex unit vec- 
tors dog = (0,0,1), d^ig = ^(l,-i,0) and d+ig = 
-^(1, i, 0). We assume a real polarization basis, e^^ 0. 

In the Rotating Wave Approximation only the terms 
that conserve the total number of excitations in the 
atomic ensemble and the quantized field are retained in 
the interaction Hamiltonian Eq. ([SJ 



RWA 



N 

i=i k,\ 



ik-ri 



(gj . \ at e"*'"''"-' 
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The RWA Hamiltonian Eq. ([7]) couples only the states 
connected by the horizontal arrow in Fig. [1] Therefore, 



if there is initially only one excitation in the system, the 
wave function at any later time has the form 

N 1 

+ E^£aW«"^^^*I5)|1£a)- (8) 

k,X 

where we have introduced the following notations; \g) = 
ffi: ff2, 5Ar) with no atomic excitations, and |ep = 
gi, ej, ■■■,gN) with one excited atom, respectively. |0) 
and |ljr_^) denote photon states with no photon and with 
one photon in the k mode with the polarization A, re- 
spectively. 

Substitution of Eq. ^ into the time depen- 
dent Schrodinger equation with the Hamiltonian from 
Eqs. dH), (g]) and © yields 

N 1 

= E E (4. • %0 (9) 



3 = 1 i'=-l 



and 



A" - - E es,3.e-((— )-^^--) {d,g . . (10) 



In the Appendix |X] we derive explicitly all the matrix 
elements needed for the analysis. 

Starting with an atomic excitation, and thus the ini- 
tial condition e^„{0) = 0, the Equations ©-([TU]) can be 
formally integrated. Here we additionally assume the 
Markovian approximation, i.e., there is no atomic popu- 
lation dynamics faster than the time scale of the atomic 
decay: /3J(r) « /3J(t) for |t - r| < E"! [H. This yields 

N 1 * 

(11) 

Substituting the last equation into Eq. ([TU| . we obtain a 
closed set of equations for (3^: 



N 1 



V 2 ik-(ri-rj) , 



/^.''^-EEE/^j^^e 
j=i ^^=^1 k,\ 

X {d,g ■ ej,) (4. . 6^,) dre-'(— )(^-*) (12) 

Since the polarizations of the photons enter into 
Eq. (jl2l) as separate scalar products, we directly perform 
the sum over the polarizations using fll 

E {"^ng ■ Hx) (4^ • %a) = '^vg ' - fc^) " dg,.- (13) 
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After introducing the short hand notations Rij = n— -fj , 
C-qu = drjg ' ( ~ '^gi^j the resulting set of equations 
for the atomic excited state amphtudes becomes 



- Jo 



fc 



i^" k 

^^g-H'^o-<^.)(r-t)^ (14) 

where we have exphcitly separated the terms with j — 
I and V — rj since these terms have different physical 
interpretations. 

The next step in our derivation involves evaluation of 
the sum over field propagation directions. To this we 
pass to the free space continuum using the standard form 
Y,k ^ /o°° '^^k'^l I di^{k) [9*1, and we arrive at 



V 



(27rc) 



Jo Jo 



X / dijJkUjlgle'^-^'^ / dTe-'('^°-"'=)(^-*) ] . (15) 
Jo Jo 

The angular integration can now be 

performed noting that ^'^^/C^) 



2^ 

47r| / and 



Efce''=-^/(fc), Jdn{k)[^I ~kk 

J dn{k)e~''''^ (^f - kk^ = AirVikR) with k = 

Rij = \Rij\ and the second rank tensor ^{kR) defined 
as [13 



ikR) = ( 



kR 



(16) 



This results in a set of equations 

" (27rc)3n 
2 



+ ^/3J(l-5,,<5,,)x 
X / dujkUjlD^^, f dTe-*('^«-"'=)("-*^ ] (17) 



with Dl^i = d^g ■ "^{kRlj) ■ dg^. 

We now turn to the integral over time using the stan- 
dard expression [l^ [l^ 

where P denotes the principal value. 

The set of Eq. (fT7|) are now in the form of linear cou- 
pled equations for the excitation amplitudes of the dif- 
ferent atoms 



r 



A 



E E (19) 

3 = li^=-l 



We have explicitly grouped terms where couples to 
itself and where it couples to different /3J as they have 
different physical interpretations and consequences. 



In particular, we recover F 



4d.i„k^ 



the single atom 



decay rate, in the first term in Eq. ([191) . The second term 
in Eq. (HH) with 



A = P 



dujki^l 

OJk - i^O 



(20) 



has imaginary value and does not depend on the loca- 
tion or number of the atoms in the ensemble. This terms 
represents the single atom Lamb shift. It formally di- 
verges, also when the RWA is not applied, but since it is 
a constant, we assume that it can be dealt with by the 
conventional assumption that it is properly included in 
the measured atomic transition frequency. 

Furthermore, in the second line of Eq. we have 
separated the contributions which explicitly depend on 
the geometry of the atomic sample. The third term is 
real valued and describes the modification of the decay 
rate of the atom I much in the analogy to the single 
atom effects. Unlike the single atom decay, the third 
term stems from the coupling to all atoms via the field 
modes and, hence, describes the collective effect: 

/"OO 

Bj = (1 - SjiSjj,y) / dwkujl5{ujk - '^o)x 
Jo 

X (d^g-^ {kRlj)-dg,y (21) 

By carrying out the Wfe integration and using 



(1 - 5jAj^) (^drjg ■ ^ (kRij) ■ dgi 



= (1 - 5ji) [d,,g ■ V {kRij) ■ dg,j , (22) 

due to the orthogonality of d^g for different Zeeman tran- 
sitions, Ea. (pT|) can be simplified to the expression 



= (1 - 6ji) ( d„g ■ V (koRij) ■ dg. 



(23) 
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where ^ is given by Eq. (ITBl) and fco = wq /c- 

Finally, the last term of the second line of Eq. 
with imaginary coupling coefficients 



X P 



djjg ■ V (fci?ij ) • dg 



(24) 



describe effective level shift of the state |ej') caused by 
the dipole coupling to all the other atoms, mediated by 
the quantized radiation field. 

It has been observed in the case of scalar and vec- 
tor photons interacting with an ensemble of two level 
atoms, that the RWA Hamiltonian yields an incorrect 
description of the single atom Lamb shift and the col- 
lective coupling of the atoms [1, |^. The discrepancy 
stems from neglecting the contributions from the electric 
dipole coupling terms which are present in Eq. ©, but 
suppressed in the RWA interaction Hamiltonian Eq. ([7]) 
and which drive the off-resonant, virtual transitions de- 
picted schematically with the diagonal dashed arrow in 

Fig.m 

For the case of two-level atoms it was observed that 
these contributions can be correctly represented by re- 
taining only the RWA Hamiltonian, if one extends, rather 
artificially, the lower integration limit of ujk to " — oo" 
when we go to the free space continuum in Eq. (ITSi) Q . 
Using this rule, the coefficients in our equations of mo- 
tion for the case of multilevel atoms change. The single 
atom Lamb shift now becomes 



Ai^) = p 



duj 



(25) 



The decay rate contributions due to the collective ef- 
fects described by Eq. (|2T|) are unchanged, since is 
positive according to our definition. In contrast, the col- 
lective level shifts of state |ej') due to the presence of all 
other — 1 atoms change to 



G 



u(A) 



(1 



X P 



(jJk - Wo 



(26) 



The equations of motion Eq. (IT5)) with the collec- 
tive interaction coefficients in the form Eq. ((2T|) and 
Eg (pS)) are consistent with the result by Friedberg et al. 
[Tol] based on the direct implementation of the pairwise 
dipole-dipole interaction between atoms. Since the mere 
observation of the equivalence between the two results 
does not offer a solid physical argument that they should 
also be equivalent with the result of a full treatment of 
the interaction Hamiltonian, we will in the following sec- 
tion derive the equations of motion of the atomic excited 
state amplitudes taking into account the virtual transi- 
tions caused by the non-RWA terms in the interaction 
Hamiltonian (3). 



III. BEYOND THE ROTATING WAVE 
APPROXIMATION 

The derivation in the present section to a large extent 
follows the derivation in Sec. [Hi and we will explicitly 
highlight the additional terms that appear due to the vir- 
tual transitions to states with more than a single excita- 
tion. Due to the parity conservation dictated by Hv^ the 
only state coupled directly to |ep|0) is \Rn^m)\'^k\) '^ith 
two atoms excited and one photon present in the electro- 
magnetic field. Due to its large violation of energy con- 
servation this transition is significantly suppressed rela- 
tive to the energy conserving transition to and 
the coupling to states with even higher excitations via 
the state |e5^e^)|lj:_^) can be safely neglected. The state 
vector of the atoms and the quantized field can hence be 
expanded on the form [1, [l^ 

N 1 

i^w) = E E /3;we-^"°*iepio)+ 

j=l u=-l 
N 1 

+ E E E"r„£.w^"'^"^~^'""'*i«)ii*;A)- 



m— n+1 ' 



(27) 



where in addition to Eq. ^ we have introduced a 
state |e^ei^)|lg^) = |.gi, e^, ej;,, gAr)|lg^) with 
two atoms excited and one photon in the field mode. 

The substitution of Eq. ([27]) into the time dependent 
Schrodinger equation with the Hamiltonian from Eq. ([T]) 
and Eq. ^ yields: 



E ^kX 



gke 



i((uk~uja)t-k-ri) (j . \ 



k,\ 
N 



y y ya"^ - gk 



"=1 , ''.M=-l fe.A 
m—n+1 



Ak-Trr 



(28) 



N 1 



= E E /JJs^e'""'-'')*-^^-^^) (dg, ■ e>) (29) 

3=1 f=-l 



and 



+ /3;g,e^«-' (dpg ■ e>) (30) 



The coupling to the atomic ground state jg), 
which has exactly the same form as in the RWA 
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case, and the coupling to states with two excited 
atoms are obtained using the fact that there are 
only four nonzero matrix elements with the Hamilto- 
nian Eq. ®: (0|(e7|7?„t|5)|l,-;,), {l^^^H.^nAe'im . 
(0|(e;'|i7.„t|e-e(;)|l,^J and {l^^^\{e':^el\H,nt\eim . We 
present the explicit form of these matrix elements in the 
Appendix [X] 

Since the mathematical structure of the new terms 
in Eq. ((28)) - (p0)) is formally similar to the correspond- 
ing terms for excitation number conserving transitions, 
we integrate Eq. (pg|) -([5n | using the initial conditions 
650- (0) = and ^^^^^.(O) = 0, as in SecHH Here again we 
additionally assume the Markovian approximation, and 
obtain 



N 1 

EE 



(31) 



+I3;{t)e 



''pg ^qry 



Eq. p3l) as separate scalar products, we directly perform 
polarization summations according to Eq. (I13p . After in- 
troducing short hand notation Onm = if m < n and 
Omyi = 1 if m > n, the set of equations simplifies 



N 1 



^l'--T.T.j:^^9l f dre 



N 1 



E T.j:^l'9lf dre' 

1 .._ 1 r. "'0 



Hu:o~Uk){T-t) . 



m— 1 fi——l 



(34) 



Regrouping the terms and noting that 6ij + 9ji = 1 — Sij 
we expand Eq. ([Ml) : 



X <?,y dTe'("''+'^°)^ (32) 




These expressions are inserted into Eq. (|28p to yield a 
closed set of equations for the single-atom-excitation am- 
plitudes f3f: 

N 1 

/3? = -E E E/5j3'^*^""'^^x 

J=i "=-1 k,X 



N 1 

E E E 

m=n+l '^'^ 



^9 ^kX 



Jk-{r,„-r„)g A -L 



+ Pm [ dgv • er 



kX 



{d'^g • %a) Sim5r,p. ] ■ (33) 



The first coupling term is equivalent to the energy con- 
serving contributions in Eq. pUj) in the RWA, while the 
subsequent terms come from the coupling to states with 
different excitation numbers. Even though the popula- 
tion of these states is negligible due to their large vio- 
lation of energy conservation, they perturb the coupling 
between the states with a single atomic excitation. In 
order to interpret their physical meaning we perform a 
few mathematical manipulations of the non-RWA cou- 
pling terms. Since the polarizations of photons enter into 



/3;'--/3;'E5' c,, f 

Jo 



Cr,r, / dTe-*''^°-""''(^-*)- 



Cuu{^ - Sij) (1 - 5i,j5u,r^) / dre 



i{uo+ijJk)(j — t) 



J. I' k 



^ik-Rij / ^^g-i(wo-"fc)(-r-t)^ 



(1 - (5ij)e'^--^^' / dre-''("«+"'')(^-*) 



, (35) 



where we have explicitly separated the terms with j — I 
and i> = rj to write the equations in a form similar to 
Eq. (HID. 

The first term in the refers to single atom effects, i.e., 
the real part yields the decay and the imaginary part 
yields the Lamb shift of the atomic excited state. Com- 
pared with the RWA analysis, which yields the first co- 
efficient in the square bracket, the coupling to the two- 
excitation sector provides a correction to the single atom 
coupling term. The second term couples the excited state 
amplitudes of the different atoms to each other and is re- 
sponsible for the collective effects. Also here, the first 
term in the square bracket is obtained in the RWA, while 
the second term is the non-RWA correction to the col- 
lective coupling terms. As in the previous section, we 
can introduce the free space continuum of field modes 
as in Eq. (jl5l) and perform the angular integration as in 
Eq. ([T7]), which yields 
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Noting that from the definition of we have 



• dgi. ^ dgi, - "^ {kRij) ■ d. 



- 5ij){l - SijS^n) \ dre 



-i(Wo+'^):)(T-t) 



+ 



Jo 

^^g-i(wo-Wfc)(T-t)_|_ 



(36) 



Finally, carrying out the time integrals and using 
SiLi - ^j7)(l ~ '^ji'^r,^) = iV - 1, we arrive at 



the equation Eq. ((T9)) . but with the modified coefiicients, 



UJk - UJQ 



{N - 1)P 



duk^l 



(37) 



The first term yields the RWA contribution to the Lamb 
shift of the excited atomic state found also in Eq. (l20t . 
The second term is due to the virtual transition with 
any one of the other [N — 1) ground state atoms being 
excited, i.e., it represents the single atom ground state 
Lamb shift of those atoms 0-Q ■ 

The coupling of different atomic excited state ampli- 
tudes in Eq. (|19l) explicitly depends on the geometry of 
the sample. The real part describes the modification of 
the decay rate of the atom I due to all the other atoms 
around it 



B 



u(F) 



(1 - (o 



(fco-Rij 



(38) 



and is unchanged, relative to the RWA case. The virtual 
transitions to the higher excited states thus turn out to 
have no direct influence on the decay terms. 

On the contrary, the imaginary coupling terms in 
Eq. (|19l) . which yield dispersive coupling among the ex- 
cited state atomic amplitudes and collective level shift ef- 
fects, acquires an extra contribution, relative to the RWA 
case in Eq. (1^ . 



X P 



UJk - OJQ 



dr,g ■ ^ {kRij) ■ dg 



X P 



{dgy ■ ^ {kRij) ■ dn 



(39) 



and using Eq. (j22l) , we can simplify Eq. (1391) and combine 
the terms with the two frequency denominators 



G;(^) = (1-5,,)x 



X P 



oo - t^O 



Y{kRij)-dgA, (40) 



We have thus recovered exactly the same expression as 
one obtains from the RWA ansatz and the purely formal 
extension of the frequency integral to negative values in 
Sec. im In the Appendix 151 we evaluate the last integral 
explicitly. 



IV. CONCLUSION 

Using the Rotating Wave Approximation and the full 
Hamiltonian, we have presented parallel derivations of 
the equations of motion describing the atomic population 
dynamics in a system of N multilevel atoms interacting 
with the quantized vector electromagnetic field. We find 
agreement of the full method with the results of the RWA, 
when the latter is extended by contributions from neg- 
ative frequency photon states. A similar agreement was 
found for simpler atomic transitions and for scalar elec- 
tromagnetic fields. The RWA with the simple extension 
to negative frequencies is easier to handle in practical 
calculations, and also easier to generalize to more com- 
plicated problems, but neither the simplest cases nor our 
more general analysis, offer a clear explanation of its va- 
lidity. With reference to Fig. [51 one may advocate that 
the virtual transitions governed by the non-RWA terms 
in the full Hamiltonian explore the same frequency range 
and degree of energy violation. This is, indeed, what 
we see explicitly in passing from Eq. to Eq. (|40p . 
The same figure, however, shows that the two virtually 
excited state spaces differ by the kind of atomic states 
involved, and, in particular by the very different dimen- 
sionality accompanying the combinatorial aspects in se- 
lecting two rather than zero excited atoms. The fact that 
the off resonant coupling to the excited states in the left 
part of Fig. [2] leads to the same reduced dynamics in 
the single-excitation subspace, as the artificial coupling 
to the negative energy continuum, must be due to the 
detailed counting of coupling terms, responsible for the 
similar equivalence observed by Milonni et al. (9|]. The 
agreement may still be surprising, given the fact that the 
actual states coupled experience the excited state Zee- 
man degeneracy, while the negative frequency continuum 
describes only ground states. The rotational invariance of 
the coupling to all field propagation directions and polar- 
izations must be at the root of the apparent suppression 
of any azimuthal quantum number dependence of the re- 
sult. The latter may thus be affected if the system is 
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FIG. 2: (Color online) Quantized matter light interaction treated by two diflerent methods. The left panel shows how the full 
Hamiltonian Eq. ((Gjl, resonantly couples the states {|ej)|0)} to the states A)} and off-resonantly (dashed line) to the 

higher excited states {|eJ^e{Ji) A)}. The right panel shows the same resonant coupling and artificial ofF-resonant coupling to 
states with negative photon energy. In this letter we show that the two methods yield identical corrections to the interactions 
in the resonant manifold of states. 



for example enclosed by non-trivial boundary condition 
for the electromagnetic field, such as in the vicinity of 
dielectric surfaces or mirrors. 

From a practical perspective, one may argue, that the 
virtual transitions, which we neglect in the RWA, have 
only very small implications for the collective superradi- 
ant emission for the ensemble The collective 

decay, however, also features subradiant states and, 
e.g., states with atoms excited with a position depen- 
dent phase that favors emission in directions which are 
suppressed by the dipole selection rule of the individual 
atoms [T^]. These phenomena may be sensitive probes 
for the small effects addressed in this manuscript. 



and 

m^Mr.C\H,nt\eim = {me'i\H„^tKe>:^)\lj:^))^ ■ 

(A4) 

Appendix B 

Since ^ given by Eq. ([T6l) consists of three parts, we 
integrate each part separately using [l7l | 
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Appendix A 



and 



The matrix elements with single atomic excitations are 
evaluated straightforwardly 



{0\{e'^\H,nt\ 



-ihgkC 



ik-ri 



(d,g ■ e-jj (Al) 



and 



(1. 



(A2) 



The matrix elements with two atoms excited are evalu- 
ated after noting that t?^e|ej:^e^) = ii j n,m and oth- 
erwise we can formally write (T™|ej:^e^) = \en)rf"l\e!^): 



SimSnti) (A3) 



„ sin(A:i?) , 
P I dk . \ ' = TTCOs{kR) 



k — ko 



(B3) 



with k — (jJk/c and ko — loo/c. Consequently, defining 



V(A:i?) 



Y - rr) 

J kR 



Eq. (1301) can be rewritten as 



G) = {1- 5,1 ) (d„g ■ V {kRij ) • dg,) ■ (B5) 
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